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ABSTRACT 

Early -type galaxies (ETGs) are known to possess a number of quite useful scaling 
relations among their photometric and/or kinematical quantities. We propose a unified 
picture reducing both the fundamental plane and the photometric plane to suitable 
projections of a single relation. Modelling the ETG as a two component system, made 
out of a luminous Sersic profile and a NFW dark halo, and applying the virial theorem, 
we are able to express the velocity dispersion oo as a function of the effective intensity 
(I e ), the effective radius R e and the Sersic index n. In a log -log plot, this relation 
reduces to a hyperplane (i.e., a plane in a four dimensional configuration space) which 
we dubbed the Sersic Virial Hyperplane (SVH). The tilt of the SVH can be fully 
explained in terms of a power - law scaling of the stellar (rather than the global) 
mass - to - light ratio T* with the total luminosity Lt, while the scatter is determined 
by those on the c — M V i r relation between the concentration c and the virial mass 
M V i r of the dark halo component. In order to test whether such the observed SVH 
is consistent with the theoretical assumptions, we perform a detailed simulation of 
a large ETGs sample reproducing the same photometric properties of a SDSS low 
redshift ETGs catalog. It turns out that the simulated SVH is fully consistent with 
the observed one thus validating our theoretical approach. 

Key words: galaxies : elliptical and lenticular, cD - galaxies: kinematics and dynam- 
ics - galaxies: fundamental parameters - dark matter 



1 INTRODUCTION 

Elliptical and SO galaxies (hereafter, collectively referred to 
as early -type galaxies, ETGs) present a striking regularity 
in their luminosity distribution. The ETG surface brightness 
is well fitted by the well known r 1,/4 profile (de Vaucouleurs 
1948), while considerable better results are obtained using 
the Sersic r 1 ^" law (Sersic 1968). From the photometric 
point of view, therefore, ETG may be considered as char- 
acterized by only three parameters, namely the slope n of 
the Sersic profile, the effective radius R e containing half of 
the total luminosity, and the effective surface brightness fi e 
defined as /J, e — n(R e ), or equivalently the average surface 
intensity (J e ) = (L T /I) / (jxR\). 

The ETG kinematic may be schematically character- 
ized through its central velocity dispersion <to which, un- 
der suitable assumptions on the mass profile, gives informa- 
tion on the mass and hence the mass - to - light (hereafter 
M/L) ratio. Being a mass tracer, it is reasonable to expect 
that (To is somewhat correlated to the total luminosity Lt, 
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even if it is difficult to forecast an analytic form for such a 
correlation, given the subtleties of the luminous and dark 
components modelling. It is therefore not surprising that 
empirical searches for such a correlation were early under- 
taken. A remarkable success was represented by the Faber- 
Jackson (FJ) Lt tx <tq relation (Faber & Jackson 1976). The 
large scatter in the FJ law led to the need for higher dimen- 
sional representations of ETGs. Considering n = 4, there re- 
mains just three parameters describing an ETG so that one 
could wonder whether a relation exist among the photomet- 
ric quantities (R e , {Ie}) and the kinematical ao- This relation 
were indeed found (Djorgovski & Davis 1987; Dressier et al. 
1987) and, when expressed in a logarithmic scale, is just a 
plane soon dubbed the fundamental plane (FP). It is worth 
noting that such a plane was not unexpected. Indeed, a sim- 
ple application of the virial theorem gives R e oc Oq{I^)~ , 
under the hypotheses of constant M/L ratio and structural 
homology. The observed FP plane is however tilted, i.e. one 
indeed finds R e oc aS(R)- , but with (a,/3) = (1.49,0.75) 
rather than (2, 1) as forecasted before (Bernardi et al. 2003). 
Such a tilt may be easily explained introducing a power - law 
correlation M/L oc L T (with 7 ~ 0.27), but interpreting the 
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origin of such a relation is a difficult and ambiguous task 
due to proposals ranging from non - homology (Prugniel & 
Simien f997), to stellar populations effects (Dressier et al. 
1987), from systematic variations in kinematical structure 
(van der Marel 1991; Bender et al. 1992; Busarello et al. 
1997) to a combination of different terms (Trujillo et al. 
2004). 

Although the de Vaucouleurs profile is a satisfactory 
fit, it is well known that the Sersic profile has to be pre- 
ferred (Caon et al. 1993; Graham & Colless 1997; Prugniel 
& Simien 1997). As such, forcing n — 4 in the fit may sys- 
tematically bias the estimate of R e and (I e ) and hence affect 
the FP. Introducing n increases the number of parameters 
needed to describe ETGs leading to wonder whether scaling 
relations exist. Actually, given the observational difficulties 
in measuring ao, it is worth looking for empirical correla- 
tions involving only the photometric parameters. Interest- 
ing examples are the Kormendy relation (KR) between R e 
and He (Kormendy 1977) and the scalelength - shape rela- 
tion between R e and n (Young & Currie 1994). However, 
just as the FJ relation is a projection of the FP, both the 
R e - /i e and the R e -n relations may be seen as projections 
of a more fundamental law among these three photometric 
parameters. In logarithmic units, such a relation indeed ex- 
ists and it is a plane referred to as the photometric plane 
(PhP) recently detected in both near infrared (Khosroshahi 
et al. 2000) and optical (Graham 2002). While observation- 
ally the PhP is confirmed also at intermediate redshift (La 
Barbera et al. 2004, 2005), a definitive theoretical interpre- 
tation is still lacking. Modelling the stars in ETG as a self- 
gravitating gas, Lima Neto et al. (1999) have recovered a 
PhP like relation (referred to as the entropic plane) by as- 
suming that the specific entropy (i.e., the entropy for mass 
unit) is constant for all ETGs. Later, Marquez et al. (2001) 
derived an energy - entropy (or mass - entropy) line giving a 
possible explanation for the structural relations among pho- 
tometric parameters. Moreover, they also find out that the 
specific entropy increases as a consequence of merging pro- 
cesses so offering a possible way to test the model against 
the observed variation of the PhP with redshift. 

There are two general lessons to draw from the above 
short summary. First, two dimensional scaling relations turn 
out to be projections of a more general three dimensional 
law. It is therefore worth wondering whether this also holds 
for the FP and PhP being possible projections of a four pa- 
rameters law. It is worth noting that a first step towards 
this direction has been attempted by Graham (2002) fitting 
a hyperplane (in logarithmic units) to the set of parame- 
ters (n, R e , (I e ) , (To), but it was never prosecuted. On the 
other hand, from a theoretical point of view, both the en- 
tropic plane and the FP are tentatively explained on the 
implicit assumption that ETGs are in a state of dynamical 
equilibrium so that the virial theorem applies and the Boltz- 
mann - Gibbs entropy may be evaluated. Motivated by these 
considerations, here we investigate whether a four dimen- 
sional relation among photometric and kinematic quantities 
may come out as consequence of the virial theorem and some 
assumptions on the stellar M/L consistent with stellar pop- 
ulation synthesis models. Should such a relation exist, one 
could thus reconcile both the FP and the PhP under the 
same theoretical standard thus representing a valid tool to 
investigate ETG formation theories. 



As a first key ingredient, we describe in Sect. 2 the as- 
sumed mass models for both the luminous and the dark com- 
ponents of ETGs. In Sect. 3, we first compute the quantities 
entering the virial theorem according to the model detailed 
above and then obtain the four dimensional scaling relation 
we are looking for. In order to test whether our derivation is 
consistent with what is observed, we simulate ETG samples 
that are as similar as possible to a large SDSS based sample 
following the procedure extensively discussed in Sect. 4. The 
results of such a testing procedure are described in Sect. 5, 
while Sect. 6 is devoted to fit our theoretical relation to the 
observed data using a Bayesian approach. Discussion and 
conclusions are finally held in Sect. 7. 



2 MODELLING ELLIPTICAL GALAXIES 

Although some recent results show the presence of thin discs 
in inner regions of elliptical galaxies^ , ETGs may be well 
described as a luminous stellar distribution embedded in a 
dark matter halo dominating the outer mass profile. In the 
following, we will assume spherical symmetry for both these 
components. While this is an acceptable hypothesis for the 
halo, it is clearly an oversimplification for the elliptical lu- 
minous component. Nevertheless, this will allow us to get 
analytical expressions for the main quantities we are inter- 
ested in without dramatically affecting the main results. 



2.1 The Sersic profile 

Under the hypothesis of constant M/L, the surface density 
of the luminous component may be easily obtained as : 



E(R) = T+I(R) 



(1) 



with Y* the M/L ratio and I{R) the surface luminosity 
density. As well known, the Sersic r 1 ^™ law (Sersic 1968) is 
best suited to describe the surface brightness distribution of 
elliptical galaxies (Caon et al. 1993; Graham & Colless 1997; 
Prugniel & Simien 1997). Motivated by these evidences, we 
will therefore set : 



I(R) = I e exp i —b n 



(I) 



l/n 



(2) 



with I e the luminosity density at the effective radius R e and 
b n a constant defined such that the luminosity within R e is 
half the total luminosity. It is possible to show that b„ may 
be found by solving (Ciotti 1991) : 



F(2n,fo n ) = T(2n)/2 



(3) 



where T(a,z) is the incomplete F function and T(a) the 
actual r function. Useful approximating formulae may be 
found in Graham & Driver (2005) and references therein, 
but we will exactly solve Eq.(3) in the following. 

Assuming cylindrical symmetry, the luminosity profile 
within R is : 

t SO galaxies contain, by definition, a thin disc, so that, strictly 
speaking, our following discussion applies only to the bulge com- 
ponent. However, neglecting this disc does not introduce any sig- 
nificant systematic error. Moreover, our sample is mainly domi- 
nated by elliptical galaxies so that we confidently neglect the disc 
component in SO systems. 
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L{R) = 2tt 



f I(R')R'dR' 
Jo 



Lt x 



7(2n, b n x) 
T(2n) 



with x = R/R e and 

L T = 2-KnI e R 2 e b- 2n c bn Y{2n) 



(4) 



(5) 



the total luminosity. The volume luminosity density j(s) 
may be easily obtained deprojecting I(R). Defining s — 
r/R e (with r the radius in spherical coordinates), we get 
(Mazure & Capelato 2002) : 



v{s) = -- 



1 f 00 di dx 



dx ^ 



with i(x) = I(R)/I e and z/(s) = (R e /I e )j(s). Some algebra 
finally leads to the following expressions for the mass density 
p*(s) and the mass profile M*(s) : 



MSj ~ 4nRl I M {s) ' 
M*(«) = Ml x Tm(s) 



(6) 
(7) 

(8) 



Im(s) 

where we have denoted with Ml the total stellar mass : 
Ml = 4(6 n /n)T*7 e ^T M (oo) , 
having used the abuse of notation 
/(co) = lim f(y) . 

V— ►oo 

Finally, to get Eqs.(6) and (7), we have used the auxiliary 
functions : 

f°° x^- n) l n exp \-b n (x 1/n - 1)1 



(x 2 - S 2 ) 1/2 



(9) 



1m{s) 



[ l v (s' 
Jo 



)s' 2 ds 



(10) 



Both these functions cannot be analytically expressed, but 
are straightforward to be numerically evaluated. 



2.2 The dark halo 

Most of the kinematical tracers of the total gravitational po- 
tential are usually affected by a severe degeneracy between 
the luminous component and the dark one so that there are 
different dark halo models able to fit the same data for a 
given stellar mass profile. It is therefore important to rely 
on a physical theory of halo formation to select models which 
are both compatible with the data and also physically well 
motivated. From this point of view, numerical simulations 
of galaxy formation in hierarchical CDM scenarios are very 
helpful since they predict the initial shape of the dark mat- 
ter distribution. Here, we assume a NFW profile (Navarro 
ct al. 1997) as initial dark matter halo and neglect the effect 
of the baryons gravitational collapse. The main features of 
the NFW model are : 



Pdm{v) 



Ps 



X (1+ x)' 2 



x — r/r 



Mom (f ) = 4 W f f(x) = M vir f(x)/f{c) 
f(x) = \n(l + x) 



1 + x 



(11) 
(12) 
(13) 



C — fvir I f*s i 

M - 4nSth 3 

IVJ-vir — Q pcrit^vi 



(14) 
(15) 



where c is the concentration parameter, M v i r the virial mass 
and r V ir the virial radius^. The model is fully described by 
two independent parameters, which we assume to be c and 
M v ir- Numerical simulations, supported by observational 
data, motivate a correlation between c and M V i r so that 
the NFW model may be considered as a single parameter 
family. Following Bullock et al. (2001), we adopt : 



c = 15 - 3.3 log (M^IO 12 ^ 1 M q ) 



(16) 



with a log normal scatter 51ogc ~ 0.11. 

The NFW model is not the only model proposed to fit 
the results of numerical simulations. Some authors (Moore 
et al. 1998; Ghigna et al. 2000) have proposed models with 
a central slope steeper than the NFW one. On the other 
hand, it is also possible that the inner slope does not reach 
any asymptotic value with the logarithmic slope being a 
power - law function of the r (Navarro et al. 2004; Cardone 
et al. 2005) or that the deprojected Sersic profile also fits the 
numerical dark matter haloes (Merritt et al. 2005; Graham 
et al. 2006a; Graham et al. 2006b). However, the difference 
between all these models and the NFW one is very small for 
radii larger than 0.5%- 1% the virial radius so that we will 
not consider models other than the NFW one. 



3 THE VIRIAL THEOREM 

Elliptical galaxies are known to be characterized by scaling 
relations among their kinematic and structural parameters. 
In an attempt to investigate whether such empirical laws 
may be recovered under a single theoretical scheme, we can 
rely on the hypothesis of statistical equilibrium. In such an 
assumption, the virial theorem holds : 



2K + W = 



(17) 



with K and W the total kinetic and potential energy re- 
spectively. Both these quantities may be evaluated given the 
assumed spherical symmetry as we detail in the following. 



3.1 Kinetic energy 

Neglecting the net rotation velocity of the system (which 
is reasonable, given the low values of v c in elliptical galax- 
ies), the total kinetic energy is given as (Binney & Tremaine 
1987) : 



K = 3-k E(R)a 2 p (R)RdR 



(18) 



with £(.R) the star surface density and u P (R) the luminosity 
weighted velocity dispersion projected along the line of sight. 
For a spherically symmetric system, assuming an isotropic 
velocity dispersion tensor, it is : 



t The virial radius is defined such that the mean density within 
r v i r is 8th times the critical density p cr it- According to the 
concordance ACDM model, we assume a flat universe with 
(Q m , CIa, h) = (0.3, 0.7, 0.72) where 8 th = 337. 
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a 2 (R) 



?*I(R) J R 



p*(r)GMtot(r)^r 2 — R 2 



dr , (19) 



with Mtot(r) = M*(r) + Mdm(t) the total mass profile. As 
a first step, it is convenient to split the integral in Eq.(19) 
in two terms separating the dark halo contribution from the 
luminous one. It is then easy to get : 



As a next step, we insert Eq.(21) into Eq.(18) and then split 
the integral in two terms, the first one originated by the 
product E(s) x T„(s), and the second one due to S(s) x 
T° M (s). By using Eq.(22), we finally obtain: 



f 

JR 



f 

J R 



p4r)GM4r)Vr 2 - R 2 



dr = 



Ml 



l v (s)l M {s) 



K = -Mjffg x fc(p) 

with p denoting the set of parameters 

p = (n, R e , I e ,T*,M vir ,c) , 

where we have set : 



(24) 



/' 

J R 



s 2 (s 2 - x 2 y 1/2 

M?M vir 
4ti- j R|Xm(co)/(c) 

I v (s)f(s,R e /r s ) 
r , s 2 (s 2 - x 2 )~ 1/2 



ds 



p*(r)GM DM (r)Vr 2 - R 2 



dr = 



Hp) - ^ 



fc*(n) = 



k*(n,Re,I e ) + Mv J^* k° M (n,R e /r s 

3 

T£(x, n)xdx , 



(25) 
(26) 



ds 



where f(s,R e /r s ) is obtained by replacing x with (R e /r s )s 
into Eq.(13). Adding the two terms above and using Eqs.(2) 
and (8) , we can finally write : 

2 26„GMj 
°^ R) - (20) 
Tl(x,n) + (M vir /M?) f- 1 (c)I° M (x,n,R e /r s ) 
cxp [-&„ (a; 1 /™ - 1)] 

with x = R/ R e and : 

T„(s)1m(s) 



Tl{x,n) 



s 2 (s 2 — x 2 ) 



ds 



T° M (x,n, Re/r a ) 



I4s)f(s,Re/r s ) 
s 2 (s 2 — x 2 )^ 1 ^ 2 



ds . 



As we will see later, <Jp(R) does not enter the applications 
we are interested in. It is rather its central value u% averaged 
within a circular aperture of radius R e /8 which is measured 
by the galaxy spectrum. We therefore evaluate : 



(To 



/ p*(r)M tot (r) [T {R e /8r) - F (Q)] dr 
2/ flc/8 p*(r) [A (Re/8r) - A o (0)] rdr 
with (Dalai & Keeton 2003) : 

y\/l — y 2 + arcsiny — tt/2 y<l 



(21) 



My) 



A (y) 



y > i 



arcsin y y < 1 
7r/2 y > 1 



When evaluating Eq.(21), we make the simplifying assump- 
tions Mtot — M+ since the dark halo mass is negligible with 
respect to the stellar one over the range interested by the 
integral. By such an assumption, it is easy to get : 

gmT 



<T = 



Re 



x lo(n) 



having defined : 



I (n) 



1 j^l v {s)I M (s) [^ (l/8a) - ^o(O)] d t 
J ™ (°°) 2 f*' /8 T v (a) {Ac ( l/8a) - Ao (0)] sds 



(22) 



.(23) 



lo(n)lM (co, jQ 

k? M (n,R e /r s ) = ±— / l? M (x, Re/r s )xdx .(27) 
' 2o(n)iM(oo) J 

Note that, although I e and Y* do not explicitly enter the 
above equations, they are however included as parameters 
since they determine the total stellar mass Mj because of 
Eq.(8). On the other hand, r s is not counted as an indepen- 
dent parameter since it is determined as function of M„; r 
and c. Should we use Eq.(16), the virial mass M vir would 
be the only parameter related to the dark halo properties. 

3.2 The gravitational energy 

The computation of W may be carried out in a similar way 
starting from the definition (Binney & Tremaine 1987) : 



'/ 

Jo 



W = ~^G / ptot(r)M to t(r)rdr 



(28) 



Splitting the total density and mass as the sum of the lumi- 
nous and dark components, after some algebra, one gets : 



W = -2— X W(p) 



(29) 



with p denoting the same set of parameters as above. The 
dimensionless quantity w(p) is defined as: 

w(p) = to*(n) 

' M v i r \ ( R, 



+ 



(M vir y ( Re \ , . 

( M v i r \ DM i n i \ 

\M T ) W * ( c ' n '- Re / rs ) 
("^r) w *DM(c,n, R e /r s ) . 



(30) 



It is then only a matter of algebra to demonstrate that 
to*(n) 



^mM Jo 



wdm(c) = 



/ 2 (c) Jo 



lv(s)lM(s)sds , 

In (1 + cy) - cy(l + cy)^ 1 
(1 + cy) 2 



(31) 



dy , (32) 



DM, r, I N 

w+ (c,n,Re/r s ) 



1 



X M (co)/(c) 

/ l v (s)f(s,R e /r s )sds , (33) 



poo 

/ M 

Jo 
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■w* DM (c,n,R e /r s 



Xm(oo)/(c) 
R, 



J q (i + j L sJl M (s)sds . (34) 



It is worth noting that, while n and R e directly enter the 
integrals above, T* and 7 e only work as scaling parameters 
through MT . Moreover, a qualitative analysis shows that 
the first term in Eq.(31) turns out to be the dominating one 
so that w(p) is essentially a function of n only. 

3.3 Scaling relations from the virial theorem 

Inserting Eqs.(24) and (29) into Eq.(17) and solving with 
respect to ao, we get : 

2 GMj w(p) 

^ Re fc(p) ' 

It is then convenient to introduce : 
L T /2 



■kRI 



(35) 



so that the total stellar mass is Mj = T*Lt = 2-KT+(I e )R 2 e 
and the relation above can be recast as : 



2 log do = log (I e ) + log R e + log 



Hp) 



+ log (2ttGT,) . (36) 



Let us suppose for a while that it is possible to neglect 
the dark halo component. Should this be the case, both 
fc(p) and w(p) turn out to be a function of (n,R e ,{I e )), 
where hereafter we use (I e ) rather than I e as a parameter^. 
Although determining how the ratio tu(p)/fc(p) depends 
on (n,R e ,(Ie)) needs for a full computation of the inte- 
grals involved, we can, as first approximation, suppose that 
w(p)/k(p) is linear in a logarithmic scale. It is possible, 
therefore, to write : 

log Hp)/fc(p)] ~ a log {Q + b log Re + clog (n/4) + d .(37) 

From stellar population synthesis models, we know that the 
stellar M/L ratio may be correlated with the total stellar 
luminosity Lt- Approximating this relation as a power - law, 
we can therefore write : 

log T* ~ a + 13 log L T = a + (3 log (2n{I e )R 2 e ) . (38) 
Inserting Eqs.(37) and (38) into Eq.(36), one finally gets: 
log do = a T log (Ie) + b T log R e + c T log (n/4) + d T (39) 
with 

a + /3 + l 



aT = 



6t 



CT 



b + 2/3 + 1 



(40) 



dT = | + ^log(2^) + ilogG + d 



3 See, e.g. Graham & Driver (2005), for the relation between Ie 
and (Ie) and other related formulae. 



As a first remark, let us note that the exact value of d de- 
pends on the adopted units. In the following, we will express 
(Jo in km/s, (I e ) in L /pc 2 and R e in kpc. In particular, hav- 
ing expressed R e in linear rather than angular units makes 
d dependent on the distance to the galaxy. A second impor- 
tant caveat is related to our starting hypothesis of having 
neglected the dark halo component. Actually, we do know 
that galaxies are embedded in their dark matter haloes. As a 
consequence, w(p)/fc(p) is a function of the halo parameters 
too. To take into account this dependence, we still assume 
Eq.(37), but let d be an unknown function of (c, M v i r ) to be 
determined by the data. 

With all these caveats in mind, Eq.(39) defines an hy- 
perplane in the logarithmic space allowing to express the 
kinematic quantity logoo as a function of the photometric 
parameters log(/ e ), \ogR e , log(n/4). Since we have recov- 
ered it for a Sersic model under the hypothesis of virial equi- 
librium, we will call it the Sersic Virial Hyperplane. 

Should our assumptions hold for real elliptical galaxies, 
the Sersic Virial Hyperplane (hereafter SVH) should repre- 
sent a tight scaling relations among kinematic and photo- 
metric parameters. For an idealized sample of galaxies per- 
fectly satisfying our working hypotheses and all at the same 
distance, the scatter around this hyperplane should be gen- 
erated by essentially two terms. First, the halo parameters 
(c,M v i r ) may differ on a case -by -case basis. This is the 
same as saying that the dark matter content in the inner 
regions or, equivalently, the global M/L ratio (defined as 
Mtot/Lr rather than Mj/Lr)is different from one galaxy 
to another. On the other hand, dr in Eq.(40) may scatter 
from one galaxy to another because of different values of 
the parameters (a, j3) of the Y* — Lt relation because of 
different details of the stellar evolution process. Note that 
this latter effect may also affect the coefficients (ot, br) thus 
further increasing the scatter on the SVH. 

It is worth stressing that the SVH may be seen as a 
generalization of both the FP and PhP which, from this 
point of view, reduce to particular cases of the SVH. Indeed, 
forcing the de Vaucouleurs model to fit the galaxies surface 
brightness profiles is the same as setting n = 4 in Eq.(39). 
Solving with respect to log R e , we get : 



log Re = a FP log a + b FP log (Ie) + c FP 
which is indeed the FP with 
a F p = 1/br 



(41) 



b F p 



c F p 



-aT /br 
-dr /br 



(42) 



Actually, we do not expect that the coefficients of the ob- 
served FP are equal to what is predicted by Eq.(42) since, 
assuming n = 4 in the fit, can bias the estimate of (R e , {Ie)) 
in a way that depends on what the true value of n is. More- 
over, the departure of n from n = 4 introduces a further 
scatter which is not included in the above expression for 
c F p. Although these effects have to be carefully quantified, 
it is nevertheless worth stressing that the FP turns out to be 
only a projection of the SVH on the plane log (n/4) = so 
that its coefficients may be (at least, in principle) predicted 
on the basis of physical considerations. 
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On the other hand, when solving Eq.(39) with respect 
to log -R e , we can also neglect the dependence on logao as- 
suming this latter is, in a rough approximation, the same 
for all galaxies. We thus get : 

log R e = a PhP log (J e ) + bp hP log (n/4) + c PhP (43) 
with 

aphP = —ar/br 



bphP = —c T /b T 

cphP = (1/6t) log cr — dr/br 



(44) 



which is indeed the PhP. Note that, since log ao is obviously 
not the same for all galaxies, the scatter in the PhP may 
then be easily interpreted as a scatter in log ao and hence in 
the kinematic structure of the galaxies. 



4 TESTING THE SVH 

The derivation of the SVH in Eq.(39) relies on two main hy- 
pothesis which are analytically formalized in Eqs.(37) and 
(38). While the relation T* oc Lj, may be tested resorting 
to stellar population synthesis models, checking the validity 
of Eq.(37) needs for a detailed computation of w(p)/fc(p) 
as function of the photometric parameters (n, (I e ),R e ), the 
stellar M/L ratio Y* and the halo parameters (c, M v i r ). Per- 
forming such a computation over a fine grid in this six di- 
mensional space is prohibitively expensive. We can however 
rely on a different and more reliable approach. Rather than 
evaluating to(p)/fc(p), we simulate a sample of galaxies with 
given values of the above parameters and use the set of equa- 
tions in Sect. 3 to compute ao and the ratio w(p)/k(p). We 
thus end up with a sample of simulated measurements of 
(ao, n, (I e ), R e ) which we fit with the SVH thus determining 
the coefficients (a sim , b sim , c s i m , d sim ). Let us then assume 
that our simulated sample have the same distribution for the 
parameters (n, {I e ),R e , T*) as a real sample. Let us then de- 
note with (a obs , b obs , c obs , d obs ) the values obtained by fitting 
the SVH to the observed sample. Should our assumptions be 
satisfied, we have to find : 

&obs — Oj S im "I - [3/2 



bobs — b s i m -\- (3 
Cobs — Csim 

d ob s — dsim + a/2 + {(3/2) log (2vr) 



(45) 



with (3 an estimate of the slope of the T* — Lt relation. 
Moreover, the scatter in the observed SVH should be the 
same as the one evaluated by the simulated sample. Actually, 
since we do not have a model independent guess for (3, we 
can estimate (3 as 2(a o6s — a sim ) or as b obs — b S i m - The two 
estimates thus obtained should of course be equal thus giving 
a further check on the validity of the model. 



4.1 The data 

In order to carry on the approach detailed above, a key step 
is a sample of ETGs as large as possible. To this aim, we 



have started from the NYU Value - Added Galaxy Catalog 
(hereafter, VAGC) which is a cross - matched collection of 
galaxy catalogs maintained for the study of galaxy formation 
and evolution (Blanton et al. 2005) and mainly based on 
the SDSS data release 4 (Adelman - McCarthy et al. 2006). 
Among the vast amount of available data, we use the low - 
redshift (hereafter, lowZ) catalog of galaxies with estimated 
comoving distances in the range 10 < d < 150 h~ 1 Mpc. 
We refer the reader to Blanton et al. (2005) and the VACG 

website^ for details on the compilation of the catalog 1 1 . 

We shortcut the lowZ catalog only retaining those data 
we are mainly interested in and rejecting all the galaxies with 
no measurements of ao leaving us with 24387 out of 28089 
objects with magnitudes in the five SDSS filters u'g'r'i'z' . 
In order to select only ETGs, we apply a set of criteria which 
we briefly details below. 

(i) A Sersic profile has been fitted to the surface bright- 
ness profile of each galaxy using an automated pipeline re- 
trieving the parameters (n, R e , A) with R e in arcsec and A 
the total flux in nanomaggies. As a first criterium, we select 
only galaxies with 2.5 < n < 5.5, where the upper limit is 
dictated by the code limit n = 6.0. This selection is per- 
formed using the fit in the i' band since it is less affected by 
dust without the reduced efficiency of the z' band. 

(ii) As a second criterium, we impose the cut Rgo /R50 > 
2.6 (Shimasaku et al. 2001) with Rgo and _R 5 o the Petrosian 
radii containing 90% and 50% of the total luminosity as es- 
timated by the data. Note that, although the Petrosian radii 
do non depend on any fitting, the ratio Rgo /R50 is correlated 
with n so that the two criteria are somewhat redundant. Re- 
moving one of them or changing the order does not alter in 
a significant way the final sample. 

(iii) We exclude all galaxies with ao < 70 km/s since the 
dispersion velocity measurements for these systems may be 
problematic (Bernardi et al. 2005). 

(iv) Elliptical galaxies are segregated in a well defined 
region of the color - magnitude plane. In order to further re- 
strict our sample, we therefore impose the cut (g — r)_ < 
9 — r < (g — r) + with (g — r)± = pM r + q ± 8. Here, M r 
is the absolute magnitude in the r filter and the parame- 
ters (p, q, S) have been tailored from Fig. 2 in Bernardi et 
al. (2005) where a different ETG sample has been extracted 
from the SDSS DR2. 

The final sample thus obtained contains 5142 galaxies out of 
an initial catalog containing 24387 objects. It is worth noting 
that most of the rejected objects have been excluded by the 
first three cuts (retaining only 5172 entries), while the fourth 
cut only removes 30 further galaxies. This is reassuring since 
the last cut is somewhat qualitative and based on a different 
set of selection criteria** (Bernardi et al. 2005). 

We then use the data reported in the lowZ catalog for 



* http : //sdss .physics .nyu. edu/vagc/ 

II Note that the version we are using is updated only to the second 
SDSS data release (Abazajian et al. 2004) covering an effective 
survey area of 2220.9 square degrees. 

** It is worth noting that we cannot use these criteria since the 
lowZ catalog does not report the parameters which the selection 
by Bernardi et al. (2005) are based on. 
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the galaxies in the above sample to collect the quantities 
listed below. 

• Photometric quantities. While the Sersic index n and 
the effective radius Re in arcsec are directly available in 
the lowZ catalog, the average effective intensity (7 e ) is not 
present. To this aim, we first convert the total flux A (in 
nanomaggies) reported in the catalog in the apparent total 
magnitude m t as (Blanton et al. 2005) : 

ra t = 22.5 - 2.5 log A . 

We then use the assumed concordance cosmological model 
to estimate the total absolute magnitude Mt as : 



Mt = m t -51ogD L (2:) + 51ogft-101og(l + z) 
- K(z) — Aq — 42.38 

where z is the galaxy redshift, Dl the luminosity distance, 
K(z) the K - correction, Aq the galactic extinction, and the 
term 10 log (1 + z) takes into account the cosmological dim- 
ming. While K(z) and Ag are reported in the catalog for 
each of the five SDSS filters, our use of the luminosity dis- 
tance makes the estimate of Mt cosmological model depen- 
dent. However, for the values of z involved, the dependence 
on the cosmological model is actually meaningless. We fi- 
nally estimate : 



10" 6 x ^ [{Mt 



■M e )/2.5] 



2-kRI 



(46) 



with dex(:r) = 10 x , Mq the Sun absolute magnitude in the 
given 

filtertt. 

We stress that, in Eq.(46), R e is expressed in 
kpc rather than arcsec. To this aim, we simply use : 

i? e (kpc) = R e (arcsec) x Da(z)/206265 

with Da(z) the angular diameter distance in Mpc. 

• Kinematic quantities. The lowZ catalog reports the ve- 
locity dispersion and its error as determined from the SDSS 
spectrum of the galaxy. This is measured in a circular aper- 
ture of fixed radius R ap = Rsdss = 1-5 arcsec, while o~o in 
Eq.(22) has been estimated for R ap = -Re/8. To correct for 
this offset, we follow j0rgensen et al. (1995, 1996) setting: 



obs lowZ 

o~o = Oo x 



/ Rsd 



ss 



V tfe/8 



(47) 



with o" °™ the value in the catalog and R e in arcsec here. 

• Auxiliary quantities. The lowZ catalog contains a 
wealth of information on each object that we really do not 
need for our analysis. We do, however, add to our catalog 
some further quantities that we will use for check. In particu- 
lar, we include the absolute magnitude Msdss as estimated 
from the image rather than the fit, and the average effective 
surface brightness computed as (Graham & Driver 2005) : 

(Me) = (tle)abs + 10 log (1 + z) + K(z) + Aq 

with 

{tie)abs =M t + 2.5 log (2nR 2 e ) + 36.57 

tt We use Mq = (5.82, 5.44, 4.52, 4.11, 3.89) for the u'g'r'i'z' fil- 
ters respectively as evaluated from a detailed Sun model reported 
in www.ucolick.org/~cnaw/sun.html 



with R e in kpc. Note that (/u e ) rather than log (I e ) is often 
used in the FP and PhP fit. 

Although the Sersic law is known to well fit the surface 
brightness profile of ETGs, it is worth noting that our deriva- 
tion of (7 e ) relies on extrapolating the fit results well beyond 
the visible edge of the galaxy. As such, it is possible that Mt 
provides a biased estimate of the actual total absolute mag- 
nitude of the galaxy which is better represented by Msdss- 
A bias in the estimate of Mt propagates on the estimates 
of the colors which may be related to the stellar mass by 
population synthesis models. In order to reduce as more 
as possible such a bias, we have studied the histogram of 
Acol = col obs - colest, where col obs = Msdss,] ~ MsDSS.k 
and colest = Mt,j — Mt,k- In principle, all these histograms 
should be centered at the null value with a small scatter. 
After removing 42 outliers, this is indeed the case for g' — r' 
and r' — i' (with rms values of 0.08 and 0.06 mag respec- 
tively), while this is not for u' — r' and i' — z' (having rms 
values of 0.21 and 0.13 mag). Motivated by this bias, we will 
use only data in g'r'i' filters when fitting the SVH taking, 
in particular, the results from the fit to the i' band as the 
fiducial ones. 



4.2 The simulated sample 

The approach we have outlined above to test the validity of 
the SVH (and, hence, of the hypotheses it relies on) is the 
construction of simulated sample of ETG which is as similar 
as possible to the real one. In particular, the photometric 
parameters of the simulated sample should match as closely 
as possible those of the observed one. 

To this aim, we first look at the histograms of the Sersic 
parameters (n, R e " , log (I e )) with R e " the effective radius 
in arcsec, the effective radius R e in kpc, and the distance 
modulus dm computed as : 

dm = 51ogDz,(z)-51ogft+101og(l + ,z) 
+ K(z) + Aq + 42.38 . 

Note that all the quantities entering the definition above 
are available for each galaxy in the lowZ catalog so that the 
distribution of dm values is easily obtained. We use the i! 
values as fiducial ones and start from the histogram thus 
obtained to compute the cumulative distribution functions 
(CDF) for the quantities of interest. A simulated galaxy is 
then generated according to the steps we sketch below. 

i. Using the CDFs obtained above, we generate the i' val- 
ues of the parameters (n, R e " , log {Ie), Re, dm). It is worth 
stressing that the ratio R e " /Re provide an estimate of the 
angular diameter distance Da(z) to the simulated galaxy 
and, hence, of the redshift provided a cosmological model 
has been set. As an alternative approach, one could gen- 
erate the redshift z and R e to infer R e " ■ Although the two 
approaches are in principle equivalent, we have checked that 
our choice is more stable from a computational point of view. 
Note also that generating dm directly makes it possible to 
avoid deriving an estimate for K(z) and Ag thus reducing 
the number of quantities to be generated. 

ii. The parameters we have generated are given in the i! 
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filter** , while we would like to simulate the galaxy properties 
in all the observed filters. To this aim, we have verified that 
linear relations between the parameters in different filters 
exist in the real sample so that we can write : 

p f = Apv 

with pf the set of parameters (n, R e " , log (7 e ), R e , dm) in 
the / filter and A. a diagonal matrix. Using this relations, 
we then generate the values of p/ for the simulated galaxy 
in the remaining u'g'r'z' filters. We are now able to compute 
all the described photometric and auxiliary quantities in the 
case of the real sample. 

iii. As a last step, we need to compute the velocity dis- 
persion <To and its error considering the ETG modelling de- 
scribed in Sect. 2. To this aim, we need to generate three 
further quantities, namely the stellar M/L ratio T* and the 
NFW parameters (c, M v i r ). As a first step, we construct the 
CDF for T* from the real sample. To this aim, for each 
observed galaxy, we estimate (Fukugita et al. 1998) : 



4.0 + 0.38 \t{z) - 10] 



(48) 



with t(z) the age (in Gyr) of the galaxy at redshift z com- 
puted assuming a formation redshift zf = 2. We then con- 
vert to the i' band M/L as : 



logT* = logTl +0.4 



(49) 



with V -i' = 0.79 (Fukugita et al. 1995). The CDF thus 
obtained is then used as seed for the random generation of 
the simulated galaxy T*. We can then estimate the total 
stellar mass Mj as 2%T i,R^{I e ) . In order to set the halo 
model parameters, we first generate the fraction of the to- 
tal galaxy mass represented by the dark matter. Defining 
r\ — M V i r / (M^ + M v i r ), we extract r\ from a Gaussian dis- 
tribution centred on 0.90 and with dispersion 0.03. It is then 
a simple algebra to find out M V i r = tj/(1 — rj) x Mj so that 
we have only to choose a value for the concentration c. To 
this aim, we use Eq.(16) to set the centre of a Gaussian 
distribution for c with a dispersion of 0.11 dex. Since the 
model is then fully assigned, we may resort to Eq.(22) to 
estimate ao and to Eqs.(24) and (29) to compute k(p) and 
w(p). The error on the predicted ao is finally set as sao with 
e extracted from a CDF obtained by the same distribution 
as the observed one. 

We repeat the procedure outlined above N — 10000 times 
thus obtaining a simulated sample to which we apply the 
same selection procedure as for the real data. We finally 
end up with a simulated sample containing (approximately) 
the same number of ETGs as the real one. While, by con- 
struction, the distribution of the photometric parameters 
(n, R e , (J e )) is the same as the observed one in the i' filter, 
we have checked that the same applies in the other filters. 
Moreover, we have checked that the statistical distributions 
of mt and Mt and of the colors are the same as the observed 
ones. We are therefore confident that the simulated sample 
is a fair representation of the real data thus making sense to 
use it as input for testing the SVH as described above. 



tt Note that dm depends on the filter used because of the terms 
K(z) and Aq that are wavelength dependent. 



5 RESULTS 

The real and the simulated samples described above are the 
input ingredients for the procedure used to test the SVH. 
As a preliminary step, however, we have to choose a method 
to fit the SVH to a given dataset. As well known (Feigel- 
son & Babu 1996; Akritas & Bershady 1996; La Barbera et 
al. 2000), fitting a linear relation in a multiparameter space 
may seriously depend on the method adopted to get the es- 
timate of the coefficients. Moreover, the choice of the most 
suitable method also depends on the uncertainties on the 
model parameters. In the present case, there is a further 
complication since the lowZ catalog does report the mea- 
surement uncertainties on ao, but not on the photometric 
parameters. As such, we have therefore decided to consider 
ao as the dependent variable and neglect, in a first approx- 
imation, the measurement errors. We therefore use a direct 
fit approach, i.e. we minimize : 



[( lo S ff °)i, &« - Oog^o), 



SVH} 



(50) 



where (log cro)^ obs and (log <To) iest are respectively the ob- 
served value and the predicted one through the SVH (39) 
for the i - th galaxy and the sum is over the Nobs elements in 
the (observed or simulated) sample. We estimate the param- 
eters (a, b, c, d) by minimizing the x 2 defined above with an 

iterative 1.5a - clipping rejection scheme^. The root mean 
square of the fit residuals is then taken as an estimate of the 
intrinsic scatter around the SVH. 

An important remark is in order. Although one may 
question on what fitting procedure is the most suitable one, 
we stress that we are here mainly interested in testing the 
consistency between what data effectively mean and what 
our theoretical assumptions predict. To this aim, it is more 
important that both the observed and the simulated samples 
are examined in the same way rather than discussing on 
which method is the best one. We are therefore confident 
that the choice of the direct fit approach does not alter the 
main conclusions we will draw from the comparison among 
the observed and predicted SVH coefficients and scatter. 

As a final warning, note that in the rest of this sec- 
tion, we only consider data and simulations carried out in 
the i! filter in order to avoid any bias from our qualita- 
tive generation of the galaxy properties in other filters. We 
have however checked that the main results are unchanged 
should we use the g' or r' filters. Moreover, we only report 
numerical values of one fiducial simulated sample. However, 
averaging over different simulations has no significant effect 
on the results. 

5.1 The linear relation hypothesis 

As a first fundamental hypothesis in deriving the SVH 
from the virial theorem, we have assumed that the ra- 
tio w(p)/fc(p) is log-linear in the photometric parameters 
(n, R e , (I e )). With the simulated sample at hand, we may 

The tj - clipping procedure retains more than 85% of the data. 
The direct fit with no a - clipping leads to consistent estimates of 
the SVH parameters, but the scatter turns out to be artificially 
increased because of clear outliers. 
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1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 

a X] + b x 2 + c x 2 + d 



Figure 1. On the j/-axis, logcro is reported as observed for the 
ETG sample. On the x - axis, the same quantity is reported as fit- 
ted by the linear relation ax\ + bx2 + CX3 + d with xi = log (7 e ), 
X2 = logi? e , X3 = log (n/4). The coefficients (a,b,c,d) are esti- 
mated by a direct fit as explained in the text. 

check whether this is indeed the case. To this aim, we have 
used the direct fit approach fitting Eq.(37) to the values of 
log [w(p)/fe(p] computed for the galaxies in the simulated 
sample. We find a linear relation among y = log [w(p)/fe(p)] 
and x — a log (J e ) + blog R e + clog (n/4) + d thus validating 

Eq.(37). For the simulation reported, we find" : 
(a, b,c) ~ (-0.05,0.04,0.29) 

with an intrinsic scatter Omt — 0.08. Using such values and 
Eq.(40) with f3 = 0, we then estimate : 

(a T ,b T ,c T ) ~ (0.48,0.52,0.14) . 

Should our ETG modelling and the virial theorem apply, we 
then expect that the SVH in Eq.(39) with (a,b,c) fits the 
observed data reasonably well. On the one hand, we have not 
used any T* — Lt relation (that is why we have set /3 = 
above) when computing w(p)/k(p) and ao- However, stel- 
lar population synthesis (SPS) models suggest that such a 
relation indeed exists. One can, in principle, use a SPS code 
to get an estimate of (a, (3), but this should introduce a fur- 
ther set of assumptions (such as, e.g., star formation history, 
metallicity, internal absorption) needed as input for the code 
itself. In order to escape these problems, we have therefore 
decided to follow a different strategy solving Eqs.(45) with 
respect to (a,/?), using as input the observed and predicted 
values of the SVH parameters. Note that, since /3 may be 
solved from two different equations, this will offer also a 
consistency check for the theory. 

5.2 The SVH coefficients 

In order to test definitively our hypotheses, we first fit the 
SVH to the observed data using the direct fit approach out- 
lined above. The result is shown in Fig. 1 where we plot y 
vs x with y = log 00 and x as above, finding 

(a obs , b obs , c obs , d obs ) = (0.493, 0.640, 0.145, 0.754) (51) 

IH "We are not concerned now with the value of d since this also 
depends on the average distance of the galaxies sample. 




1.8 2 2.2 2.4 2.6 2.8 3 

a X! + b x 2 + c x 3 + d 



Figure 2. On the y-axis, we report logcro as computed for the 
simulated ETG sample. On the a; -axis, we instead report the 
same quantity as fitted by the linear relation ax\ + bx2 + cx$ + d 
with x\ = log(Te), X2 = log R e , X3 = log (n/4). The coefficients 
(a, b, c, d) are estimated from a direct fit as explained in the text. 

with an intrinsic scatter a°^1 = 0.045. As Fig. 1 clearly 
shows, the SVH fits the data remarkably well over the full 
observed range in log ct . The fit residuals are quite small and 
do not correlate neither with x nor with y. We can therefore 
safely conclude that the SVH is indeed observationally well 
foundedllH. 

Motivated by these encouraging result, we now turn to 
the simulated sample repeating the same fit as above, but 
with y now computed on the basis of our ETG modelling 
and assumptions on the halo parameters. Fig. 2 clearly shows 
that the SVH fits very well the simulated data with : 

{a slm ,b slm ,c sim ,d stm ) = (0.475,0.515,0.142,0.822) (52) 

with an intrinsic scatter of™ = 0.040. Comparing Eq.(51) 
and (52) is a rewarding task. Indeed, c obs and c s im are re- 
markably close which is exactly what we do expect because 
of Eq. (40) . Moreover, the scatter is quite similar with a mod- 
est underestimate (Aaint = 0"°nt — of™ = 0.05). Considering 
that we have included as source of the scatter only the one 
due to the c— M v i r relation and neglected that introduced by 
the conversion from Y]f to T*, such a small value of Acr» nt 
has not to be considered as a shortcoming of our theoretical 
assumptions. 

Solving Eq.(45) with respect to /3 using the values of 
a obB and a sim gives /3 ~ 0.03, while (3 ~ 0.12 is obtained 
from b obs and b 3 i m . These two estimates are not consistent 
with respect to each other so that one could be tempted 
to conclude that something is wrong with our hypotheses. 
There are however some important systematic errors which 
these difference could be ascribed to. First, we stress that the 
values in Eq.(52) have been obtained from a simulated ETG 



III It is worth noting that a first attempt to fit a hyperplane 
to a small sample of Coma and Fornax galaxies (collected from 
the literature available at that time) has yet been successfully 
performed in Graham (2005). However, in that work, the fit was 
only empirically motivated and the rest of the paper gives off 
the dependence on <ro concentrating on the PhP. We can there- 
fore consider the present work as the first theoretically motivated 
study of the SVH based on a large and homogenous ETG sample. 
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sample obtained under the simplifying hypothesis of spher- 
ical symmetry. As well known (Binney & Tremaine 1987), 
for a given total mass, the velocity dispersion may consider- 
ably differ between spherically and flattened systems. As a 
second simplification, we have computed the kinetic energy 
in the case of isotropy in the velocity space, while it is likely 
that ETGs are anisotropic systems. Quantifying these effects 
is not possible unless one takes carefully into account the 
distribution in the intrinsic flattening q and the anisotropy 
parameter /3 CT which are largely unknown. We can roughly 
quantify the effect of these systematics errors by first writing 
the true velocity dispersion as : 



2 

,true 



2 . 2 
°"0 T &0,sys 



al x 



1 + 



with <To, ays the term we are neglecting due to the systematic 
errors above. Using the SVH for al, we get : 

log a , true = a T log (J e ) + 6 T log 7? e + ct log (n/4) + d T 



ilog 



1 + 



/ (?0,sys \ ' 

\ ao J 



(53) 



The fit to the real data is consistent with the existence of 
the SVH so that we can argue that the last term on the r.h.s. 
is linear in logarithmic units. In particular, by writing 



log 



1 + 



/ CT0,sys \ * 

V cro / 



<51og(/ e > 



(54) 



the SVH is recovered provided that we replace or with 
ht + 8/2. The matching between simulated and observed 
coefficients is now obtained for : 



fa S im + P/2 + 8 — a bs 
bsim + /3 = bobs 

From the values of (a obs , b b s ,cisim,b S i m ), we thus get : 



(55) 



(P,5) ~ (0.12,-0.09) . 

It is worth noting that such a small value of /3 indicates a 
quite weak correlation between Y* and Lt- This is qualita- 
tively consistent with the finding of Padmanabhan et al. 
(2004). Estimating the stellar M/L from the correlation 
with the D4000 strength (Kauffman et al. 2003), these au- 
thors find T* to be almost constant with Lt which com- 
pares well with our estimated /3. Moreover, defining the 
dynamical mass Md yn as the total mass within (approxi- 
mately) R e and estimating it as Md yn = (1.65) 2 alRe/G, 
Padmanabhan et al. finds Mdyn/Lr oc L T A7 . With our find- 
ing T* oc L T 12 , we can recover the Padmanabhan et al. 
result provided Mdyn/M^ oc L^ ' 04 . Using the real data 
and estimating Y* from Eq.(48), we find M dyn /M+ oc L^ 03 
so that our model predicts Mdyn/Lr oc L^ 15 in remarkable 
agreement with the Padmanabhan et al. value. 

Notwithstanding the encouraging result for f3, it is 
worth stressing that Eq.(54) is actually an unmotivated as- 
sumption. Nevertheless, we can qualitatively assess whether 
such a hypothesis is reasonable by a qualitative analogy with 
spiral galaxies. Some recent works (Roscoe 1999; Noorder- 
meer et al. 2007) claims that the shape of the rotation curve 
is determined not only by the total luminosity, but also 
weakly depends on the mean surface density. Considering ao 



as the analog of v c for ETGs, we may qualitatively assume 
that the power -law relation — Lt describes the depen- 
dence of a ,true on Lt, while Eq.(54) takes into account a 
possible effect related to (J e ). Although this interpretation 
is, strictly speaking, unmotivated, the low 8 value and the 
positive match with the Padmanabhan et al. finding make 
us confident that Eq.(54) is not unrealistic. 

Considering these encouraging results, we therefore con- 
clude that the theoretically predicted SVH can reasonably 
coincide with the observational one. 



6 THE OBSERVED SVH 

In the previous section, we were mainly interested in testing 
the theoretical bases of SVH so that what is pressing is the 
need for a consistent and homogenous analysis of both the 
observed and the simulated data. Such a consideration has 
motivated our choice of the simplified fitting method which 
is, however, not best suited to derive the true SVH coef- 
ficients. In particular, we have not taken into account the 
measurement uncertainties on cro giving the same weight to 
each ETC in the sample. Moreover, the scatter in the SVH 
has been estimated a posteriori, while one should more cor- 
rectly take into account such an intrinsic scatter a priori, 
i.e. as a further parameter to be determined by the fitting 
procedure. 

The Bayesian probabilistic approach offers an ideal 
route to solve this problem. We do not enter in any de- 
tail here referring the interested reader to the vast literature 
available (see, e.g., D' Agostini (2004) and refs. therein). Let 
us suppose that a linear relation holds as : 



y = axi + bx2 + cx-a + d 



(56) 



and let o» n t be its intrinsic scatter. If the errors on the vari- 
ables involved are statistically independent, one can demon- 
strate that the best fit estimate of the parameters (a, b, c, d) 
and of the scatter ai n t is obtained by minimizing the follow- 
ing merit function (D Agostini 2005) : 



1 H 



AT 



1 \ ^ (Vi — axi,i — bx2.i — cx-i.i — d) 2 



■(57) 



with a Vt i and a^i the errors on y and Xj for the i- th object 
and the sum is over the M objects in the sample. It is worth 
stressing that the minimization with respect to d may be 
performed analytically, i.e., for given (a,b, c, a int ), the best 
fit d is given by : 



7 'Lt+ a l,,+ a ' 2 ' 7 'i,i+ b2 ' T 2^+ c2a 3,, 



\pAf _ 1 



(58) 



With this value of d, one can compute — In C and then find 
the set of parameters (a, b, c, d, ai„t) that minimizes it, there- 
fore representing the best fit solution. 

The general formulae (56) - (58) may be easily adapted 
to our problem. Eqs.(39) and (56) may be identified setting : 

y = log a , xi = log (J e ) , x 2 = log R e , x 3 = log (n/4) . 
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Table 1. Best fit coefficients and intrinsic scatter of the SVH in the g'r'i' filters obtained as described in the text. 
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Filter abed <Ti n t 

g' 0.433 ±0.008 0.595 ± 0.010 0.102 ±0.014 0.880 ± 0.022 0.089 ± 0.002 

r' 0.460 ±0.008 0.600 ± 0.010 0.113 ±0.014 0.841 ± 0.022 0.081 ± 0.002 

i' 0.469 ±0.008 0.594 ±0.010 0.146 ±0.014 0.824 ± 0.022 0.077 ± 0.003 



Table 2. Best fit coefficients and intrinsic scatter of the SVH as rewritten in Eq.(59) in the g'r'i' filters. 



Filter 


a e 


be 


c e 


d e 




g' 

r' 
i' 


0.897 ± 0.020 
0.973 ± 0.022 
1.020 ±0.020 


-0.659 ± 0.007 
-0.678 ± 0.008 
-0.695 ± 0.007 


0.064 ± 0.022 
0.076 ± 0.023 
0.052 ±0.021 


0.052 ± 0.039 
-0.111 ±0.039 
-0.189 ±0.036 


0.110 ±0.002 
0.100 ±0.002 
0.100 ±0.003 




Note that, for the data at hand, o\ = oi = az = 0. Finally, 
we find the results summarized in Table 1 for the SVH in 
the g'r'i' filters (excluding therefore the v! and z' data for 
the problems hinted to in Sect. 4) where the errors on the fit 
parameters have been estimated by 1000 bootstrap resam- 
pling. Fig. 3 gives a visual impression of the quality of the 
fit superimposing the best fit curve to the i' filter data. 

Comparing the results in Table 1 for the i' filter and 
those in Eq.(51), we find a good agreement for both c and d 
coefficients, while a less good but still reasonable matching 
is obtained for a and b. From this point of view, both the 
direct fit and the Bayesian approach seem to provide reli- 
able estimate of the SVH parameters. However, the scatter 
in Table 1 is almost two times larger than the one estimated 
by the direct fit method. This is an expected result. The two 
fitting procedures actually minimizes a similar merit func- 
tion and indeed, in Sect. 5, we have looked for the minimum 
of Eq.(57) forcing Oi n t — and setting to 1 the denominator 
in the second term. As discussed in D'Agostini (2005), this 
is likely to bias low the estimate of a int which is just what 
we find. Note that the scatter of the data around the best 
fit line is essentially the same as the intrinsic scatter so that 
it is fully explained by the sources generating (Tint, namely 
the variation in the halo parameters. 

Looking at the Table 1, it is not possible to definitively 
infer whether the SVH coefficients are filter dependent. On 



the one hand, the estimates for the different wavelengths 
are formally not consistent within each other. Nevertheless, 
the differences are so small that it is hard to say whether 
they are significant or due to selection effects or some other 
uncontrolled systematics. 

Being an immediate consequence of the virial theorem 
and our assumptions (37) and (38), we have up to now writ- 
ten the SVH as a relation giving the velocity dispersion as 
a function of the photometric parameters. However, we can 
solve Eq.(39) with respect to R e as: 

log R e = a e log a + be log (I e ) + c e log (n/4) + d (59) 
with : 
' ae = 1/6 

k de = — d/b 

which can be directly estimated for the different filters from 
the values in Table 1. A most reliable estimate is, however, 
obtained by fitting Eq.(59) to the data using the same pro- 
cedure as above thus obtaining also a value for the intrinsic 
scatter. The results are reported in Table 2, while Fig. 4 
shows the best fit line superimposed to the data in the fidu- 
cial i' filter. It is immediate to see that the values in Table 
2 are in strong disagreement with the theoretical expecta- 
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tions based on Eq.(60) and the values in Table 1. This is 
however by no means a shortcoming of the model, but only 
an expected consequence of problems with fitting an inverse 
relation so that we do not discuss anymore this issue. 

It is more interesting to note that the scatter in log R e 
is the same as the intrinsic one. Although this could also 
be somewhat a consequence of neglecting uncertainties on 
log_R e , it is an intriguing result. Indeed, while in the FP 
the scatter must be accounted for by both observational 
errors and intrinsic terms, here the intrinsic scatter is the 
only source of scatter. Moreover, such an intrinsic scatter is 
fully related to the variations in the dark matter parame- 
ters so that its origin is fully known and (to some extent) 
predictable provided a reliable c — M v i r relation is given. 

It is worth wondering how the inverse SVH (59) works 
as distance indicator. To this aim, we note that the scat- 
ter in log R e translates into a In 10 x a in t,e — 23% scatter 
in the distance estimates. Table 4 in Bernardi et al. (2003) 
reports a collection from literature of various FP determina- 
tions with the corresponding scatter on the distance ranging 
from 13% to 22% so that, taken at face values, increasing 
the number of parameters with respect to the FP (4 instead 
of 3 for the SVH vs the FP) does not ameliorate the per- 
formances as distance indicator. Some remarks are however 
in order. First, different determinations rely on different fit- 
ting methods. Should we have used the direct fit approach, 
the scatter in the distance should be reduced to a ~ 16% 
which is smaller than the 20% value obtained by Bernardi 
et al. (2003) for a SDSS sample comprising 9000 ETGs with 
0.01 < z < 0.3. A meaningful comparison should therefore 
rely on the same fitting method and we advocate the use 
of our Bayesian approach in order to not underestimate the 
intrinsic scatter. As a second issue, one has also to consider 
that our sample is made out of all ETGs in the lowZ cata- 
log without any separation in field or clusters objects. It is 
therefore worth reconsidering the fit and hence the scatter 
by performing the fit on smaller subsamples separated ac- 
cording to their environment. This is outside the aims of this 
introductory investigation and will be presented elsewhere. 



7 DISCUSSION AND CONCLUSIONS 

Early type galaxies may be considered as a homogenous class 
of objects from many point of views. A further support to 
this idea is represented by the existence of several interesting 
scaling relations among their photometric and/or kinematic 
parameters, the most famous ones being the FP (between 
log (To, log (7 e ) and logi? e ) and the PhP (where logoo is re- 
placed by the Sersic index n). In an attempt to look for 
a unified description of both these relations, we have pre- 
sented here the Sersic Virial Hyperplane (SVH) expressing a 
kinematical quantity, namely the velocity dispersion ao as a 
function of the Sersic photometric parameters (n, ,R e , {le}). 
In the usual logarithmic units, such a relation reduces to a 
hyperplane (i.e., a plane in four dimensions) where all ETG 
lay with a small thickness as inferred from the low intrinsic 
scatter a int . Just as the FJ relation is a projection of the 
FP and the KR a projection of the PhP, thus we find that 
both the FP and the PhP are projections of the SVH so that 
the scatter in these well known relations can be ascribed to 
neglect one of the SVH variables. 



Our derivation of the SVH relies on very few assump- 
tions. First, ETGs are postulated to be in dynamical equi- 
librium so that the virial theorem applies. Given that ETGs 
are old systems likely to have formed most of their stellar 
content and settled their main structural properties at z ~ 2 
(Merlin & Chiosi 2006), this hypothesis seems to be well 
founded, at least as a first well motivated approximation. 
Starting from this premise, the SVH comes out as a conse- 
quence of the virial theorem and of Eqs.(37) and (38). The 
first one relies on an approximated log -linear relation be- 
tween the model dependent quantity w(p) /k(p) and the Ser- 
sic photometric parameters, while the second one assumes 
the existence of a power - law relation between the stellar 
M/L ratio Y* and the total luminosity Lt- Given mass 
models for the luminous and dark components of a typical 
ETG, Eq.(37) has been verified by realistic simulations re- 
producing the distribution in the photometric parameters of 
a large ETG sample selected from the low redshift version of 
the VAGC catalog based on SDSS DR2. On the other hand, 
the relation oc L T is expected on the basis of stellar pop- 
ulations synthesis models and is also invoked to explain the 
FP tilt. It is worth stressing, however, that here the tilt of 
the SVH is due to the dependence of the stellar M/L on Lt, 
while, in the case of the FP, such a relation involve the global 
(stellar plus dark) M/L ratio. As a consequence, one has to 
resort to a mechanism coupling the dark and luminous mass 
density profiles in order to have M dyn / Lt oc L t , while here 
we only rely on what is predicted by stellar population syn- 
thesis models. From this point of view, therefore, the SVH 
relies only on known physics without the need of any unex- 
plained interaction between baryons and CDM particles. 

The discussion in Sect. 5 have confidently demonstrated 
that the observed SVH is consistent with our theoretical 
predictions so that it is worth wondering how the present 
introductory work can be ameliorated. To this aim, it is 
worth reconsidering the derivation of the observed SVH in 
Sect. 6. Here we have used the Bayesian approach to infer 
the estimate of the SVH coefficients and its intrinsic scatter 
from a dataset based on the lowZ catalog. Such a sample is 
however affected by known problems. First, the automated 
code used by Blanton et al. (2005) to estimate the Sersic 
parameters does not work very well. The recovered values of 
(n, R e , A) and hence of log (7 e ) are biased in a complicated 
way depending on the values of the parameters themselves. 
Modeling this bias and taking care of it in the fitting pro- 
cedure is quite difficult, but the problem should be coped 
with by checking how it affects the estimate of the SVH co- 
efficients. It is worth stressing, however, that this problem 
does not impact anyway our testing of the SVH. However, 
to test the SVH, one has to verify that, for a given distri- 
bution of photometric parameters based on a whichever real 
sample, the simulated sample and the observed one predict 
the same SVH coefficients which is indeed what we find. 
Note also that this method allows to automatically include 
the selection effects in the simulation without the need for 
accurately modelling them. The data we have used present, 
however, a second shortcoming, i.e. there is no estimate of 
the uncertainties on the Sersic photometric parameters. As 
shown in Eq.(57), these quantities enter the definition of the 
merit function to be minimized in the Bayesian determina- 
tion of the SVH coefficients. Although one can estimate what 
the effect is by artificially attaching measurement errors to 
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the input quantities, it is more convenient to solve the prob- 
lem by resorting to a different sample in order to address 
the problem hinted above. The recently released Millenium 
Galaxy Catalog (Liske et al. 2003; Driver et al. 2006) con- 
tains a detailed bulge/disc decomposition of ~ 10000 nearby 
galaxies with detailed fitting of the Sersic law to the surface 
brightness profile (Allen et al. 2006) . The code is tested and 
shown not to be biased and the errors on the photometric 
parameters are available. Cross - matching with the SDSS 
and selecting only the ETGs should provide an ideal sample 
to test the SVH retrieving a more reliable estimate of its 
coefficients and scatter thus allowing to reconsider it as a 
distance indicator. 

From a theoretical point of view, it is worth reconsid- 
ering our basic assumptions. As discussed above, the tilt of 
the SVH with respect to the virial theorem predictions may 
be fully ascribed to a power - law relation between the stel- 
lar M/L ratio T* and the total luminosity Lt- Although 
estimating the slope /3 from the tilt is, in principle, possi- 
ble, confronting with an expected value is welcome. To this 
aim, one can resort to stellar population synthesis models 
(Fioc & Rocca - Volmerange 1996; Bruzual & Chariot 2003; 
Le Borgne et al. 2005; Maraston 2005) by varying the dif- 
ferent ingredients entering the codes and looking (by trial 
and error) for the combination giving the slope /3 needed to 
reproduce the correct SVH tilt. Should these stellar mod- 
els be able to reproduce the observed colors of ETC, our 
derivation of the SVH could be further strengthened. 

A fundamental role in the ETC modeling has been 
played by the choice of the dark halo model. Although the 
NFW mass density profile is the standard one, it is never- 
theless well known that it encounters serious difficulties in 
explaining the inner rotation curves of low surface brightness 
galaxies (see, e.g., de Blok 2005 and refs. therein). Moreover, 
some recent evidences from the planetary nebulae dynamics 
have put into question the need for a significative amount of 
dark matter in the inner regions of elliptical galaxies (Ro- 
manowski et al. 2003; Napolitano et al. 2005). Estimating 
the dark matter content for the real galaxies is difficult, but 
we can rely on what we know from the modeling of the simu- 
lated galaxies since they reproduce the same SVH as the ob- 
served ones. Defining the dynamical mass* * * Md yn as the 
total mass within R, i.e. Md yn (R) = T*L(i?) + Mdm(R), 
we find a median value MDM(Re)/Md yn {Re) ^ 14% (with a 
rms value of 20%) so that inner regions turn out to be baryon 
dominated as expected. This also makes us confident that 
changing the halo model does not affect significantly the 
theoretical values of the SVH coefficients. A subtle effect, 
worth to be investigated, is the relation between the con- 
centration c and the virial mass M V i r of the halo. Actually, 
the scatter of this relation concurs in determining the scatter 
in the SVH. Said in another way, for given M v i r and stellar 
mass parameters, the scatter in the c — M vir relation intro- 
duces a scatter on MoM{Re) / Md yn {Re) thus contributing to 
the total SVH scatter. Numerical N-body simulations and 



* * * Note that here we compute the dynamical mass from the 
known values of the model parameters rather than estimating it 
from the velocity dispersion. For this reason, although the defini- 
tion is the same, our values for Md y „ differ from those one would 
obtain using the formula in Padmanabhan et al. (2005). 



semi - analytical galaxy formation models predict different 
halo models with its own c — M V i r relation and scatter. It 
should be tempting to investigate whether a large ETC sam- 
ple (free of the problems described above) could be used to 
discriminate among these different possibilities on the basis 
of the scatter they predict for the SVH. 

The aim of the present paper was mainly to introduce 
the SVH as a unifying scenario for the ETC scaling rela- 
tions. If the encouraging results presented here will be fur- 
ther confirmed, both observationally and theoretically, we 
are confident that the SVH could represent a valid tool to 
investigate the ETC properties under a single picture. 
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